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Abstract 

We study the spreading kinetics of a monolayer of hard-core particles on a semi- 
infinite, chemically heterogeneous solid substrate, one side of which is coupled to a 
particle reservoir. The substrate is modeled as a square lattice containing two types 
of sites - ordinary ones and special, chemically active sites placed at random positions 
with mean concentration a. These special sites temporarily immobilize particles of 
the monolayer which then serve as impenetrable obstacles for the other particles. In 
terms of a mean- field- type theory, we show that the mean displacement Xo(t) of the 



monolayer edge grows with time t as Xo^t) = \J ID at ln(4Z)Q,t/7ra^) , (a being the 
lattice spacing). This time dependence is confirmed by numerical simulations; is 
obtained numerically for a wide range of values of the parameter a and trapping times 
of the chemically active sites. We also study numerically the behavior of a stationary 
particle current in finite samples. The question of the influence of attractive particle- 
particle interactions on the spreading kinetics is also addressed. 

Key Words: Monolayer spreading, chemically heterogeneous substrates, dynamic per- 
colation. 



1 Introduction 



The stability and spreading kinetics of ultrathin wetting films on solid substrates are of 
technological and scientific importance in many applications ranging from coatings, paints, 
dielectric layers, thin film lubrication, microelectronic devices, to fundamental studies of 
adsorption and particle dynamics In the case of homogeneous, chemically pure 

substrates, the properties of such films are relatively well understood through a series of 
experimental and theoretical works 

However, most of the naturally occurring surfaces used in thin film experiments are 
chemically heterogeneous on nanometer to micrometer scales, e.g. due to contamination, 
cavities, uneven oxide layer, etc. On the other hand, deliberately tailored chemically 
heterogeneous substrates are also increasingly being used for engineering of desired nano- 
and micropatterns in thin films (see, e.g. Refs. [§,0). In addition, some recent studies have 
revealed a possibility of controlling the growth of biological systems by attaching them to 
structured surfaces |jTD|] and to recognize biological molecules, (e.g., proteins), selectively 
by bringing them into contact with nanostructured surfaces [|1^ . 

A considerable amount of recent theoretical, numerical and experimental work has been 
devoted to the analysis of equilibrium properties of thin films on chemically heterogeneous 
substrates. These studies focused mostly on such issues as stability of films, pattern for- 
mation, appearance of self-organized structures, as well as the impact of chemical disorder 
on the contact angle and appropriate generalization of the Young's equation [|1^-[19|. Much 
less is known, however, on spreading kinetics of ultrathin liquid films on chemically disor- 
dered surfaces. Here, the only available studies concern Molecular Dynamics simulations 
pil|,PT|] and experimental analysis |2^,|2^ of precursor films spreading on substrates with 
chemically impure sites. To the best of our knowledge, no theoretical analysis has been as 
yet performed. 

In the present paper, motivated by recent experimental studies of precursor films spread- 
ing on chemically disordered substrates |2^,|2^, we analyse the spreading kinetics of molec- 
ularly thin films on substrates with randomly placed chemically active sites. We focus here 
on systems with the so-called planar geometry, i.e. on systems in which film's thickness 
(or concentration of particles in the film in case of monolayers) varies effectively only along 
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one spatial coordinate. This typical experimental situation occurs when a solid, which may 
be a plane or a cylindrical fiber, is immersed in a liquid bath. Here, the particle concen- 
tration in the liquid film, which extracts from the macroscopic meniscus and climbs along 
the solid, varies only with the altitude above the edge of the macroscopic meniscus and is 
independent of the perpendicular, horizontal coordinate. The meniscus then serves as a 
reservoir of particles, which is in equilibrium with the spreading monolayer and "feeds" it. 
The solid substrate is modeled here in a usual fashion as a regular, square lattice of adsorp- 
tion sites; chemical heterogeneity is introduced by adding some concentration a of special, 
chemically active sites, which temporarily trap moving particles which then become ob- 
stacles for others. We analyse here, both analytically and numerically, the behavior of the 
mean displacement Xoit) of the monolayer edge. In terms of a mean- field-type approach, 
wc find that Xo(t) grows with time t as Xo(t) = y^ZD^^tlii(4D^^t7^ra^, (a being the lattice 
spacing). This time dependence, which contains a non-trivial logarithmic factor, is con- 
firmed by the numerical simulations. As well, we obtain Da numerically for a wide range 
of values of chemical sites' concentration a and of the trapping times. We also consider the 
situation when our substrate is of a finite extent along the X-axis and study numerically 
the behavior of the stationary particle current. The question of the infiuence of attractive 
particle-particle interactions on spreading kinetics is also addressed. 

The paper is structured as follows: In Section 2 we formulate a microscopic stochastic 
model of spreading kinetics. In Section 3 we derive basic equations and present their 
mean-field-type solution appropriate for situations with annealed spatial distribution of 
the chemically active sites. In Section 4 we describe our Monte Carlo simulations model. 
Results of Monte Carlo simulations of spreading kinetics of monolayers composed of hard- 
core particles and analysis of the behavior of the particle current in finite samples are 
presented in Section 5. Next, in Section 6 we consider spreading behavior in the case when 
the monolayer particles experience short-range, nearest-neighbor attractive interactions. 
Finally, in Section 7 we conclude with a summary and discussion of our results. 
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2 The Model 



As we have already remarked, our model is relevant to the following experimental situation. 
Suppose that a vertical solid wall is immersed in a bath of liquid. The liquid interface, which 
is initially horizontal, changes its shape in the vicinity of the solid wall and a macroscopic 
meniscus builds up. The size of the macroscopic meniscus (both horizontally and vertically) 
is comparable to the capillary length. After a suitable transient period an ultrathin liquid 
film (a monolayer) exudes from the static macroscopic meniscus and climbs up the solid wall 
0. In Ref. a microscopic stochastic model describing spreading kinetics of molecular 
films on chemically homogeneous, ideal substrates has been developed. Here we extend 
this approach on the situation when chemical disorder is present. 

Particles dynamics on the solid surface is generally regarded as an activated random 
hopping motion, constrained by hard-core interactions, between the local minima of a 
wafer-like array of potential wells. Such wells occur because the monolayer's particles 
experience short-range forces exerted by the atoms of the solid. Consequently, the interwell 
distance a is related to the spacing between the atoms of the substrate. Without going 
into details of the particle-substrate interactions, we suppose that for the transition to 
one of the neighboring potential wells a particle has to overcome a potential barrier. This 
barrier does not create a preferential hopping direction, but results in a finite time interval 
r between the consecutive hops, defined through the Arrhenius formula. 

To specify the positions of the wells, we introduce a pair of perpendicular coordinate 
axes {X, Y), where X is a vertical coordinate, which measures the altitude of a given well 
above the meniscus (a reservoir), while Y defines the horizontal position of this well. For 
simplicity, we suppose that the lattice of potential wells is a regular square lattice of spacing 
a (see Fig. 1). It will be made clear below that the effects we observe do not drastically 
depend on the precise form of the underlying lattice. 

Further on, we assume that the substrate contains some concentration a of immobile, 
chemically active sites, placed at random positions. For simplicity, we suppose that the 
spatial distribution of these sites is commensurate with the underlying lattice of potential 
wells, such that the chemically active sites can be viewed as occupying random positions 
on the sites of the square lattice exactly (see. Fig. 1). 
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We turn next to the definition of tlie flopping probabiiities. We suppose that the fatter 
are symmetric regardfess of whether a particfe occupies an ordinary or a chemicaffy active 
site, fn the former case, a particfe chooses a jump direction with the same probabifity equaf 
to 1/4, which means that being on ordinary site a particfe afways attempts to perform a 
hop. On the contrary, in the fatter case, there is a probabifity that the particfe stays at 
the site - a pausing probabifity e, which mirrors the chemicaf specificity of sites and hence, 
resufts in a temporaf trapping effect, fiere, the particfe sefects the jump direction with 
probabifity (1 — e)/4, where the parameter e can be expressed as 



e = l-exp(Utr/kBT'^, {!] 



ksT being the temperature measured in the units of the Boftzmann constant A;^, whife 
Utr denotes the trapping energy, Utr < 0. Note that the typicaf time r* spent by a given 
particfe being trapped by a chemicaffy active site is just r* = r/(l — e). 

Consequentfy, the site-dependent jump direction probabiiities Y) can be written 
down as 



p{X,Y) 



1/4, if the site (X, F ) is an ordinary site, 

(1 — e)/4, if the site (X, Y) is a chemicaffy active site. 



After the jump direction is chosen, the particfe attempts to hop onto the target site. The 
jump is fuffiffed if the target site is empty at this moment of time; otherwise, the particfe 
remains at its position. 

Finaffy, we view the fiquid bath as a reservoir of particfes (of an infinite capacity) which 
maintains a constant concentration Cq of fluid particfes at the edge of the macroscopic 
meniscus, i.e. the fine X = in Fig. 1 (see Ref. for more detaifs). Here, for simpficity, 
we tak;e Co = 1. The behavior for arbitrary Co wiff be considered efsewhere |^ . 



We hasten to remark; that dynamics in disordered fattice gas-type modefs, refevant to 
the one empfoyed here, has been extensivefy studied within different contexts, incfuding, 
for instance, charge carrier transport in dynamic percofating systems tracer diffusion 



within the first fayers of sofid surfaces p6| and in adsorbed monofayers tracer and 



coffective diffusion on sofid surfaces [^, in pure and disordered crystafs ^,|3TI] or 



coffective diffusion in zeofites p2|-p^. The systems anafyzed in these worlcs differ, however. 



considerably from the situation under study; here we present a first, to the best of our 
knowledge, lattice gas-type description of spreading dynamics of monolayers on substrates 
with chemical disorder. 

To close this section it might be instructive to discuss the limitations of such a non- 
interacting lattice-gas-type model. In the "real world" systems, the particles appearing on 
top of a solid substrate - adsorbed particles, experience two types of interactions: namely, 
interactions with the atoms of the underlying solid - the solid-particle (SP) interactions, 
and mutual interactions with each other - the particle-particle (PP) interactions. The SP 
interactions are characterized by a repulsion at short scales, and an attraction at longer 
distances. The repulsion keeps the adsorbed particles some distance apart of the solid, while 
attraction favors adsorption and hinders particles desorption as well as migration along 
the solid surface. In this regard, our model corresponds to the regime of the so-called 
intermediate localized adsorption ll|,|3^: the particles forming a monolayer are neither 
completely fixed in the potential wells created by the SP interactions, nor completely 
mobile. This means, the potential wells are rather deep with respect to the particles 
desorption (desorption barrier Ud ^ ksT), so that only an adsorbed monolayer can exist, 
but have a much lower energy barrier Vi against the lateral movement across the surface, 
Ud ^ Vi > ksT. In this regime, any monolayer particle spends a considerable part of its 
time at the bottom of a potential well and jumps sometimes, solely due to the thermal 
activation, from one potential minimum to another in its neighborhood; after the jump is 
performed, the particle dissipates all its energy to the host solid. Thus, on a macroscopic 
time scale the particles do not possess any velocity. The time r separating two successive 
jump events, is just the typical time a given particle spends in a given well vibrating around 
its minimum; as we have already remarked, r is related to the temperature, the barrier 
for the lateral motion and the frequency of the solid atoms' vibrations by the Arrhenius 
formula. 

We emphasize that such a type of random motion is essentially different from the 
standard hydrodynamic picture of particles random motion in the two-dimensional "bulk" 
liquid phase, e.g. in free-standing liquid films, in which case there is a velocity distribution 
and spatially random motion results from the PP scattering. In this case, the dynamics 
can be only approximately considered as an activated hopping of particles, confined to 
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some effective cells by the potential field of their neighbors, along a lattice-like structure of 
such cells (see, e.g. Refs. ||3^,|3^). In contrast to the dynamical model to be studied here, 
standard two-dimensional hydrodynamics presumes that the particles do not interact with 
the underlying solid. In realistic systems, of course, both the particle-particle scattering 
and scattering by the potential wells due to the interactions with the host solid, (as well as 
the corresponding dissipation channels), are important 38]. In particular, it has been 
shown that addition of dissipation to the host solid removes the infrared divergencies in 
the dynamic density correlation functions and thus makes the transport coefficients finite 
On the other hand, homogeneous adsorbed monolayers may only exist in systems in 



which the attractive part of the PP interaction potential is essentially weaker than that 
describing interactions with the solid; otherwise, such monolayers become unstable and 
dewet spontaneously from the solid surface. As a matter of fact, for stable homogeneous 
monolayers, the PP interactions are at least ten times weaker that the interactions with 
the solid atoms . 

Consequently, the standard hydrodynamic picture of particles dynamics is inappropri- 
ate under the defined above physical conditions. Contrary to that, any adsorbed particle 
moves due to random hopping events, activated by chaotic vibrations of the solid atoms, 
along the local minima of an array of potential wells, created due to the interactions with 
the solid [|l],^]. As we have already remarked, in the physical conditions under which such 
a dynamics takes place, the PP interactions are much weaker than the SP interactions 
and hence do not perturb significantly the regular array of potential wells due to the SP 
interactions. In our model, we discard completely the attractive part of the PP interac- 
tion potential and take into account only the repulsive one, which is approximated by an 
abrupt, hard-core-type potential. 

The question of the monolayer spreading in the case when some short-range attractive 
particle-particle interactions are present will be briefly addressed in Section 6. 
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3 Basic equations and a mean-field-type solution 



Let pt(X,Y) denote the local density of the monolayer particles at time moment t at the 
site {X,Y). This local density obeys the following balance equation 

r^^^^ = -P{X,Y) p,iX,Y) Y: (l-p.(X',n) + 

{X',Y') 

+ (l-p(X,r)) J2 piX',Y')p,iX',Y% (2) 

{X',Y') 

where (X', Y') denotes a nearest-neighboring to (X, Y) site, while the summation symbol 
with the subscript {X', Y') means that the summation extends over all nearest to {X, Y) 
sites. Note that the factors ^1 — Pt{X', Y')^ and ^1 — p{X, Y)^ on the right-hand-side of 
Eq. account for the steric constraints due to hard-core interactions and represent the 
(decoupled) probabilities that the target sites are unoccupied at time moment t. 

Equation (H) holds for all particles except for the rightmost particles for each fixed 
Y, since for the latter, by definition, the hops away of the monolayer (i.e. such that 
increase their X position to X + a) are not constrained by the hard-core interactions. Let 
now Xo{Y,t) denote the X-position of the rightmost particle in the column with fixed Y. 
Evidently, one has for pt{X = Xo(Y,t),Y) the following equation 

^ dp,{Xo{Yt),Y) ^ _p^Xo{Y,t),Y) p,{Xo{Y,t),Y) {l - Pt{Xo{Y,t),Y')) ' 

Y'=Y±a 

- p{Xo{Y,t),Y) pt{Xo{Y,t),Y) (^l-p,(Xo(r,t)-a,F)) - 

- p{Xo{Y,t),Y) pt{Xo{Y,t),Y) + 

+ (l-p(Xo(y,t),F))[ Yl PiXoiY,t),Y')ptiXoiY,t),Y') + 

Y'=Y±a 

+ p{Xo{Y,t) + a,Y) pt{Xo{Y,t) + aX)] + p{Xo{Y,t) - a,Y), (3) 

which thus has a different structure compared to Eq. (|^). Note that the last term on the 
right-hand-side of Eq. (4) is not multiplied by neither the occupation factor p((Xo(F, t), Y) 
nor by the steric factor (1 — pt{XQ{Y,t),Y)). This happens, namely, because the last 
term describes the event in which the rightmost particle, present, by definition, at the site 
Xo(t) — a, (i.e. pt{Xo(Y,t)—a,Y) = 1), hops at the vacant site Xo(t), (i.e. pt{Xo(Y,t),Y) = 
0). 
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We turn next to the mean-field-type picture assuming first that chemically active sites 
are uniformly spread along the substrate with mean density a, and p{X, Y) is a position- 
independent constant 

p{X, F) ^ ^ (4) 

An estimate of pa will be presented below. 

Then, we note that the dependences of pt{X, Y) on the X and the Y coordinates have 
quite different origins. There is a reservoir of particles, which maintains fixed occupation 
of all sites at X = 0. Consequently, we may expect a regular X-dependence of pt{X,Y). 
In contrast, the F-dependence may be only noise; the uniform boundary at the X = 
insures that there is no regular dependence on the Y coordinate and, in absence of disorder 
in the jump direction probabilities, only the particle dynamics may cause fluctuations in 
Pt{X,Y) along the F-axis. Hence, following Ref. we will disregard these fluctuations 
and suppose that the local density varies along the X-axis only, i.e. pt{X,Y) = pt{X). 
Consequently, we will have an effectively one-dimensional problem in which the presence 
of the F-direction will be accounted only through the particles' dynamics. We note finally 
that assumption of such a type is, in fact, quite consistent with experimental observations 
0, which show that in case of sufficiently smooth substrates and liquids with low volatility 
the width of the film's front is very narrow. 

Then, in neglect of the fluctuations along the F-axis the variable pt{X) can be viewed as 
a local time-dependent variable describing occupation of the site X in a stochastic process 
in which hard-core particles perform hopping motion (with a time interval r* between the 
consecutive hops) on a one-dimensional lattice of spacing a connected, at the site X = 0, 
to a particle reservoir which maintains constant occupation of this site. 

For t ^ T, characteristics of such a process are then described by the following non- 
linear system of coupled equations. The mean displacement of the rightmost particle (the 
monolayer edge) obeys: 

c?Xo(t) apo 
^ dt ~ ~ 
where pt{X) is determined by 



Pi(X = Xo(t)), (5) 



dptjX) ^ d^PtjX) 

^ dt 4 9X2 ' y'^) 
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which holds for < X < Xoit) and is to be solved subject to two boundary conditions: 

PtiX = 0) = 1, (7) 

and 



dpt{X = Xo{t)) a^p^ dpt{X) 
a T — — 



dt 4 dX 



-.MX = Xo(t))^ 

(X=Xo(i)) 



These two boundary conditions mimic, first, the presence of a particle reservoir, and second, 
show that for the rightmost particles of the monolayer the jumps away of the monolayer 
are not constrained by hard-core interactions. 

We note now that Eqs. (^) to (|^) constitute a classical mathematical problem of 
solving a partial differential equation with one of the boundaries being imposed in the 
moving frame, which is akin to the so-called Stefan problem. Its solution can be found in 
a standard way by observing that the density profiles pt{X) written in terms of the scaling 
variable uj = X/Xo(t) become stationary. In the limit t ^ r, the mean displacement of 
the monolayer edge thus follows 



Xo{t) = ^2DMn[^), (9) 
where is given by 

D„ = (10) 
In a similar fashion, one finds that the total number M(t) of particles, 

POO 

M{t) = / dXptiX), (11) 



emerged on the substrate up to time t, obeys 



M{t)r^j^, (12) 



which implies that the mean density in the monolayer slowly decreases with time 

M{t) 



Pt{X) 



Xo{t) 



^ Trln (4Dj/na^ 



(13) 
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Note that dependence of pt{X) on disorder, which enters only through the effective diffusion 
coefficient is logarithmically weak. Note also that the mean displacement Xo(t) of the 
monolayer edge grows at a faster rate than the conventionally expected pure diffusive \/t- 
law due to an additional factor ^J\\i{t); consequently, fitting of experimental curves or 
numerical results with a pure ^/t-law is meaningless since the effective diffusion coefficient 
will appear to depend on time of observation. In Fig. 2 we present numerical evidence 
of this additional logarithmic factor. In Fig. 3 we depict numerical results describing the 
behavior of Da- Analytical estimates of Da will be presented elsewhere |^ . 



We finally remark that within the employed mean-field dynamical approach, we can 
also obtain an average stationary particles current < Jpart >■ Solving Eq. (|]) subject 
to the reservoir boundary condition in Eq. (|^), as well as imposing a trapping boundary 
condition at the right edge of the substrate, pt{X = N) = 0, we find that 

< Jpart >= (14) 

i.e. the current has a Fickian dependence on the substrate's length. 

The effective diffusion coefficient Dj can be estimated within a mean-field-type approx- 
imation as follows: in the stationary state it matters actually how much time, on average, 
a given particle spends on a given lattice site. Such an average time is, evidently, 

< r >= r X (1 - a) + r* X a, (15) 

where the first term represents a contribution of ordinary sites, while the second one gives 
an average time spent by a given particle on chemically active sites. Consequently, the 
effective diffusion coefficient Dj can be estimated as 

^ 1 — e , , 

Dj = = 7 r 16 

4<r> 4rl-e(l-a) ^ ^ 

This result is, of course, exact for a = and a = 1, i.e. for chemically homogeneous 
substrates. It appears that it describes reasonably well (see Fig. 4) the numerical data for 
a ~ 1 and arbitrary e, as well for small values of e and arbitrary a. 
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4 Numerical simulations 



In our simulation algorithm, we follow closely the model defined in Section 2. We consider 
a square lattice A with linear sizes and Ly and with every lattice site {X, Y) we associate 
an occupation variable n(^x,Y) which may assume only two values {+1, 0}. The value +1 
signifies that the site (X, Y) is occupied, while means that this site is vacant. 

The initial configuration is an empty lattice except for the zeroth raw {X = 0). The 
left edge of the system X = is coupled to a particle reservoir which keeps the zeroth raw 
always occupied (Cq = 1), i.e. {n(o,y) = l.Y = l,...,Ly}. In the F-direction periodic 
boundary conditions are imposed to reduce the finite-size effects. The right edge of the 
system is coupled to an empty reservoir, so that the raw {X — L^+l} is always empty. Note 
that such a formulation allows us to study both dynamic and static characteristics. While 
studying spreading dynamics, we take sufficiently large and take care that displacement 
of the rightmost particle in each column Y is less than L^ + l. When studying the behavior 
of the stationary particle current, we focus on L^; not that large and let the system evolve 
until the density profiles in the system attain a stationary state. 

The following time-saving procedure has been implemented. At every non-normalized 
time step i a particle in the system is chosen at random. Let the particle's coordinates be 
denoted by {X,Y). Then the particle may either stay at the site {X,Y) with probability 
e(X, Y){— e, 0}, or with an equal probability, p{X, y) = (1 — e{X, l^))/4, may attempt to 
jump onto one of the neighboring sites, chosen at random. The jump is actually fulfilled if 
the target site is empty. Otherwise, the particle remains at the site (X, Y) . If the initial 
site is in the zeroth raw and if after the update the particle moves it is immediately filled 
by a particle from the reservoir and the number of particles Ni in the system is increased. 
If the initial site is in the last raw X = and if after the update the particle moves to 
X' — Lx + 1 the number of particles in the system is decreased. The time is renormalized 
according to 

ti+i^ti + -^, (17) 

where A^j+i is the total number of particles in the system at the non-normalized time (^ + 1). 
We use the averaged renormalized time in our studies of the time-dependent quantities. 
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Most of the simulations are performed for a system of size 100 x 25, 200 x 50 and 
100 X 100 in units of the lattice constant a. Larger system sizes are also considered in few 
cases. The results are usually averaged over A^^ = 2, 5, 10 different substrates and for each 
substrate A^^ = 5, 10 different runs are performed. Typical Monte Carlo simulation lasted 
1.6 ^ 2 X 10^ MCS per site. 

5 Simulation Results 

After passing through a transient regime the system reaches a stationary non-equilibrium 
state characterized by a stationary average particle current Jpart flowing through the system 
and a constant average density gradient. 

We studied here how do both, the spreading diffusion coefficient Da and the diffusion 
coefficient in the stationary state Dj, depend on the pausing probability e and on the 
concentration a of the chemically active sites. The spreading diffusion coefficient 
was determined from the time dependence of the average interface position Xo(t) before 
particles start leaving the right edge of the system. It appears that the law in Eq. (|) 
describes very well the time behavior of the average interface position not only for a = 

0, e = 0, but also for practically the whole interval of values of a and e, except at e = 1.0, 

1. e. infinitely deep trapping sites (see Fig. 2). 

For determination of the diffusion coefficient Dj in the stationary state we use Fick's 
law, Jpart = —DjVp, by measuring the average particle current, Jpart (per site), and the 
average density gradient, Vp(~ const.), in the stationary state at given pausing probability 
e and concentration a. The obtained results for the spreading diffusion coefficient are 
presented in Fig. 3 and the corresponding results for the diffusion coefficient Dj in the 
stationary state are given in Fig. 4. Curiously enough, the values found for are always 
lower then those obtained for Dj. 

We turn now to the special case when the pausing probability on chemically active 
sites is e = 1. The specific feature of this case is that the particle, once arriving at any 
chemically active site stays there forever, serving then as impenetrable obstacle for the 
other particles. It means that in this case one has an induced percolative behavior. The 
time behavior of the average interface position for a > 0.1 is no longer fitted well by 
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the function in Eq. (j^) (one expects that here a logarithmic time behavior should take 
place) and the above mentioned method cannot be employed to determine the spreading 
diffusion coefficient Da. For given a the averaged density distribution in the stationary 
state is still constant and Vp ~ —(1 — a)/Lx. In order to get a reliable estimates for the 
studied quantities (e.g. the particle current Jpart) the demand on the computing time as 
the concentration a ^ increases significantly since longer time runs are necessary as 
well as averaging over more substrates is needed and finally also bigger systems should be 
simulated. The approximate value found for the concentration ~ 0.4 ± 0.01 at which 
the particle current Jpart (respectively Dj) turns to zero is consistent with 1 — p^, where 
Pc = 0.592746 [^] is the critical probabihty for site percolation in the square lattice (see 
Fig. 5). 



6 Monolayer of interacting particles. 

We turn finally to the case when the monolayer particles experience short-range (nearest- 
neighbor) attractive interactions. Let us consider the simplest possible case when the 
corresponding Hamiltonian is 

H = -U n{X,Y)n{X'X) , (18) 

{X',Y') 

where U {U > 0) is the constant describing the attraction between two diffusing particles 
and the summation symbol with the subscript "(X',Y')" means that summation extends 
over the sites {X',Y'), neighbouring to the site {X,Y). We still assume the activation 
mechanism for the hopping motion of the monolayer particles; that is, the probability for 
jump depends on the trapping energy of the site (X, Y) through: 

'uux,Yy 



Pjump{-^ J ^) — exp 



We take into consideration the interaction between the diffusing particles by assuming 
that the particle "feels" the other particles when choosing the direction for the jump, i.e.: 

P,„,((X, Y), (X', Y')) = 1 exp ( '^^^^^^y^'^^^^ ) ' ^ = E ^^-((^' (^'' ^')) 



{X',Y') 
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where 

AH{{X, Y), {X', Y')) = H{X', Y') - H{X, Y) 

and H{X, Y) = —U n{X, Y) J2{x' y) ^^(^'5 is the interaction energy of the particle at 
the site (X,F). 

For high enough temperatures one may, as a first approximation, try to determine the 
diffusion coefficients and Dj in the same way as it was done for the non-interacting 
system. The temperature dependence of the diffusion coefficients determined in this way 
is shown in Fig. 6. As could be seen taking into consideration the interaction between the 
diffusing particles leads to a decrease of the diffusion coefficients. For higher temperatures 
the effect is less pronounced. 

For lower temperatures another method for determining Dj should be employed. While 
for lower temperatures the time behavior of the average interface position Xo(t) is still 
reasonably well described by Eq. (|]), the density distribution along the spreading direction 
in the stationary state is no longer linear. In Fig. 7 the corresponding density distributions 
(for homogeneous substrate, a = 0, Utr = 0) are shown for three different temperatures 
for the interacting system in the stationary state. One can see that at ksT = 0.5U there 
is clearly a phase separation though there is a stationary particle current flowing through 
the system. The interface between the two phases is approximately at X = Lx/2. At 
higher temperatures, e.g. ksT = 2.5U the density distribution is getting closer to a linear 
distribution (as in the non-interacting case) but is still not linear. This system is very 
similar to the driven diffusive system introduced by Katz, Lebowitz and Spohn where 
there is a stationary particle current flowing in the system due to a bias in the transition 
rates. 



7 Conclusions 

To conclude, we have studied the spreading kinetics of a monolayer of hard-core particles 
on a semi-infinite, chemically heterogeneous solid substrate, one side of which is attached to 
a reservoir of particles. The substrate is modeled as a square lattice containing two types 
of sites - ordinary ones and special, chemically actives sites placed at random positions 
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with mean concentration a. These special sites temporarily immobilize the particles of the 
monolayer which then serve as impenetrable obstacles for the other particles. In terms of a 
mean-field- type theory, we have shown that the mean displacement Xo{t) of the monolayer 
edge grows with time t as Xo(t) = a/ 2Dat \n{ADat / vra^) , (a being the lattice spacing). This 
nontrivial time dependence is confirmed by the numerical simulations. For a broad range of 
values of a and of the trapping times of the chemically active sites (pausing probabilities) 
Da has been obtained from extensive Monte Carlo simulations. In addition, we have 
studied numerically the behavior of the stationary particle current in finite samples. We 
have observed that, curiously enough, the diffusion coefficient deduced from the analysis 
of the data on the spreading kinetics, and the one obtained from the analysis of the data 
on the stationary particle currents, Dj, are different from each other and obey Da < Dj. 
Besides, we have found that the system displays a percolation-type behavior when e = 1 
and Q! ^ Q!c ~ 0.4 ih 0.01. In this limiting case both D^ and Dj vanish. The question 
of the influence of attractive particle-particle interactions on spreading kinetics has been 
also addressed. We have observed that taking into consideration attractive interactions 
between the diffusing particles leads to a decrease of the diffusion coefficients. For higher 
temperatures the effect becomes less pronounced, as it should. Finally, we have found that 
for sufficiently strong attractions the density distribution along the spreading direction in 
the stationary state is no longer linear and that there is clearly a phase separation, though 
the stationary particle current does not vanish. 
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Figure captions 



• Fig. 1. Schematic representation of the monolayer in contact with a particle reservoir 
on a chemically heterogeneous substrate. Gray squares denote chemically active sites. 

• Fig. 2. Plot of Xo{t)/Vi versus ■\/ln{t) - numerical evidence for time-dependent 
logarithmic corrections to the mean displacement of the monolayer edge. Circles 
denote the time moment when the rightmost particles of the monolayer reach the 
right edge of the substrate, such that the finite-size effects come into play. 

• Fig. 3. The dependence of the spreading diffusion coefficient Da is shown: (a) 
as a function of the pausing probability e at different fixed concentrations a of the 
chemically active sites, sohd circles-solid hne - a = 0.1, solid up-triangles-solid line 
- a — 0.3, solid squares-solid line - a — 0.5, solid diamonds-solid line - a — 0.7, 
solid down-triangles-solid line - a — 0.9; (b) as a function of the concentration a 
at different fixed pausing probabilities e, solid circles-solid line - e = 0.1, solid up- 
triangles-solid line - e = 0.3, solid squares-solid line - e = 0.5, solid diamonds-solid 
line - e = 0.7, solid down-triangles-solid line - e = 0.9. 

• Fig. 4. The dependence of the diffusion coefficient in the stationary state Dj is 
shown: ( function of the pausing probability e at different fixed concentrations 
a of the chemically active sites, solid circles-solid line - a = 0.1, solid up-triangles- 
solid line - a = 0.3, solid squares-solid line - a = 0.5, solid diamonds-solid line - a — 
0.7, solid down-triangles-solid line - a — 0.9; (b) as a function of the concentration 
a at different fixed pausing probabilities e, solid circles-solid line - e = 0.1, solid up- 
triangles-solid line - e — 0.3, sohd squares-solid line - e — 0.5, sohd diamonds-sohd 
hne - e — 0.7, solid down-triangles-solid hne - e — 0.9. The dotted hues are the 
corresponding analytical curves given by Eq. 17. 

• Fig. 5. Percolation threshold. The plot of Dj versus o; for e = 1. Linear extrapo- 
lation of the numerical data gives the critical value oi a — ac at which the current 
vanishes equal to ccc ~ 0.4 ± 0.01. 
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• Fig. 6. The temperature dependence of the spreading diffusion coefficient Da and 
the diffusion coefficient Dj in the stationary state. Concentration of the chemically 
active sites is a = 0.5 of the trapping sites and their trapping energy Utr is taken 
equal to Utr = —0.7. For the non-interacting system {U = 0) the results for Dj 
and for are given by solid squares-solid lines and by open squares-dotted line, 
respectively. For the weakly (compared to the trapping energy) interacting system 
{U — 0.1) - solid up-triangles-sohd hue and open up-triangles-dotted hue depict, 
respectively, the behavior of Dj and D^. For {U — 0.3) - sohd circles-solid line 
define Dj, while open circles-dashed line determine the corresponding behavior of 
the spreading diffusion coefficient Dq. 

• Fig. 7. The average density distributions along the spreading direction X in homo- 
geneous systems {a — and Utr — 0) are shown for the interacting {U — 1) system 
100 X 25 (in units of the lattice constant) in the stationary state at three different 
temperatures: ksT — 2.5U sohd hue, ksT — lU dashed line, ksT — 0.5U dotted 
line. 
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Figure 2: Pesheva and Oshanin 
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Figure 3: Pesheva and Oshanin 
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Figure 4: Pesheva and Oshanin 
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Figure 5: Pesheva and Oshanin 
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Figure 6: Pesheva and Oshanin 
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Figure 7: Pesheva and Oshanin 
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